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FOREWORD 


Transverse  oscillations  of  a  relativistic  electron  beam  propagating  in  a 
beam-induced  ion  channel  has  been  theoretically  investigated  and  experimentally 
observed.  Good  agreement  between  experimental  observations  and  theoretical 
calculations  strongly  suggests  that  these  oscillations  may  be  caused  by  the  ion 
hose  (ion  resonance)  instability. 


CARL  W.  LARSON,  Head 
Radiation  Division 


We  have  investigated  the  transverse  oscillations  of  an  electron  beam 
propagating  in  a  beam-induced  plasma  channel.  From  the  comparison  of  both 
theoretical  and  experimental  results,  it  is  strongly  suggested  that  the 
transverse  oscillations  are  caused  by  the  electrostatic  transverse  coupling 
between  the  electron  beam  and  the  ion  channel,  which  is  known  as  the  ion  hose 
(ion  resonance)  instability.  *  However,  in  previous  studies,1’^  the  usual 
assumption  was  that  the  charge  neutralization  fraction  f  (channel  ions  to  beam 
electrons  line  density  ratio)  was  a  constant  from  beam  head  to  tail.  In  this 
context,  it  is  necessary  to  develop  a  new  theoretical  model,  which  can  describe 
some  important  physical  phenomena,  for  the  case  where  f  varies  from  head  to  tail 
and  the  beam  is  self-pinched  without  an  external  magnetic  field. 

In  this  report,  we  consider  an  intense  relativistic  electron  beam,  with 
Gaussian  profile  and  radius  R,  propagating  through  an  initially  neutral  gas. 

A  A 

The  axial  velocity  of  the  beam  is  Vez  =  cez.  The  beam  space  charge  is  partially 
neutralized  by  the  ions  in  the  plasma  channel  generated  by  the  beam  itself,  with 
the  channel  electrons  assumed  to  be  instantaneously  expelled  in  the  absence  of 
external  applied  axial  magnetic  field.  Secondary  ionization  is  neglected  in  the 
low  pressure  regime  considered  here.  The  neutralization  fraction,  f(t),  as  a 
function  of  time  can  he  approximately  given  by 
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for  f  <  1.  Here  and  Nb  are  line  densities  of  ions  and  beam  electrons, 

respectively,  cr  is  the  ionization  cross  section,  nn  is  the  neutral  gas 

density  (~  pressure,  P),  and  tQ  =  (ocnn)"*.  Strictly  speaking.  Equation  (1)  is 

a  crude  estimation,  since  in  reality,  the  beam  assumes  a  trumpet  profile,  i.e., 

R  is  a  rapidly  changing  function  of  time  until  the  equilibrium  condition-^ 

2  2 

f  =  1/y  +  ZTj/unc  is  achieved.  Here,  y  is  the  beam  relativistic 
factor,  T  is  the  beam  transverse  temperature,  and  v  =  2.82  x  10“^  N^fcm-1)  is 
the  well  known  Budker's  parameter.  However,  this  approximation  makes  the 
problem  analytically  tractable,  and  is  actually  quite  good  within  the  framework 
of  the  "rigid"  beam  and  channel  treatment  considered  here. 

In  the  rigid  approximation,  the  beam  and  ion  channel  are  treated  as  rigid 
but  flexible  rods.  The  transverse  displacements  of  the  beam  and  channel 
centroids,  Db  and  ,  can  be  obtained  from 


Nb*»e 


■Fb1(Db-D1>- 


ri  ui 

Nbmi  cTF  f  clF”  =  Fib(Di  ■  V’ 


for  f  <  1.  Here  m^  and  me  are  the  ion  and  electron  rest  masses,  and  F j ( x )  is 
the  electrostatic  force  exerted  on  column  i  by  column  j,  with  relative  centroid 
displacement  distance  x.  For  two  columns  with  Gaussian  profiles  and  equal 
radi i .  we  can  obtain 


NSWC  MP  86-426 


2e.e.N  N  22 

Fij(x)  =  -  ■  j"-’J  (1  -  e  x  /2R  ) 


*  e.e.N.N.  for  x  <<  R. 

»  J  »  J 


Making  use  of  Equations  (1)  and  (3),  Equation  (?)  can  be  rewritten  as 


(1F  *  c  •§2)2  °b  ■  -  V • 


tif '  ir  =  TO(0b  '  °f)' 


2  2  2  2 

where  a  =  e  N. /yn  R  =  c  IK(kA)/17yR  »  and  n  =  YHa/nK .  Noting  that  n  «  1. 

D  0  D  0  7 


and  that  z  enters  Equation  (4)  only  through  3/az,  we  seek  solution  to 


Equation  (4)  of  the  form 


D  =  D(t )  e 


-ik  z 
z 


A  p  p 

with  3D/ 3t  <<  k  c.  It  is  furthermore  required  that  k_K  <  1,  which  assures 


the  validity  of  the  rigid  treatment.  Then  by  substituting  Equation  (S)  into 


Equation  (4a),  we  obtain 


4r  -  “h 


■r"i’ 
0  ' 


which  can  be  integrated  to  give 
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Di 


+  i 


^  •'  -M 


Di  (t ' )  dt'. 


(6) 


In  Equation  (6),  6(t)  is  determined  by  s(t)  =  (at2/2tQ  -  k2c2t)/2kzc. 

A 

Substituting  Db  in  Equation  (6)  into  Equation  (4b)  gives 


3_ 

3t 


(e 


i  6(t ) 


_3_ 

3t 


-  —  t 

4.  nil  t. 


3t  ' 


=  1 


no*zc  is(t) 
?  e 


(7) 


The  full  behavior  of  Equation  (7)  is  complicated  and  can  only  be  calculated 
numerically.  The  detailed  results  will  be  presented  elsewhere.  However, 
Equation  (7)  can  be  approximated  to 

330,  32D.  nakc  * 

_J_  +  U>  __L  *  i  d.,  (8) 

3tJ  c 


without  loss  of  important  physical  contents.  Here  6'  is  defined  by  6'  = 

(at/tg  -  k^c2)/2kzc. 

For  the  region  6'  =  0,  Equation  (8)  can  be  solved  to  gi  ve  ~  e  ,  where 


i  nc^7c 

-  j  ( 1  -  i  ft)  (-/- 


1/3 


(9) 


Away  from  this  region.  Equation  (8)  can  be  approximated  by 


32D.  nok2c  D. 


(10) 
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which  indicates  stable  oscillations  for  6'  <  0,  and  purely  growing  for  6‘  >  0. 
From  this  theoretical  analysis,  we  conclude  that,  for  an  electron  beam 


propagating  through  an  initially  neutral  gas,  the  ion  resonance  mode  for  a  given 

2  2. 


axial  wavelength  kz  becomes  unstable  at  time  t^  =  kzc  t^/ct,  and  that  growing 
oscillations  given  by  Equation  (9)  can  be  observed  only  for  very  low  pressure 


2  2 

(large  t0)  regime.  In  reality,  for  f  =  t /t n  <  1/y  +  2"!^/^  ,  the  beam  radius 


R  is  large  and  rapidly  changing,  and  thus  long  wavelength  modes  such  that 


2  2  2  2 

kzc"/a  <  1/y'  +  2T^/\mc  are  stabilized  by  wall  effects  and  betatron  detuning. 


The  unstable  mode  observed  are  expected  to  be  the  longest  axial  wavelength  that 
satisfies  the  condition 


2T 


•>?«*  i 

— ~  t + — *?• 

y  \jmc 


(ID 


and  the  onset  time  of  this  instability  is 


*d  = 


2  2 

kzc  *0 


a 


(12) 


The  experiment  is  performed  with  the  Transbeam  accelerator,  which  generates 
a  nominal  700  kV,  100  kA,  100  ns  pulse  into  a  matched  load.  The  matched 
electron  beam  diode  consists  of  a  7.5  cm  diameter  planar  carbon  cathode  and  13 
micron  titanium  anode  foil  with  a  15  mm  anode-cathode  gap.  In  order  to  obtain  a 
low  current,  high  quality  beam,  a  graphite  beam  stop  with  a  2  cm  diameter  hole 


on  axis  is  placed  immediately  downstream  of  the  anode  foil.  This  allows 
approximately  4  kA  of  the  beam  to  be  injected  into  the  gas -filled  drift 
region.  Typical  diode  current  and  voltage  traces  for  this  configuration  along 
with  further  details  of  the  accelerator  design  and  operation  may  be  found  in 
Reference  4.  Also  in  Reference  4  we  have  characterized  the  beam  propagation 
in  the  ion  focused  regime  ( IFR ) .  The  radial  profile  and  root -mean -square  (rms) 
emittance  of  the  beam  have  been  measured  with  a  simple  emittance  meter  and  a 
radiachromic  film  detector.  The  beam  is  characterized  by  a  transverse 
temperature  of  35  keV  and  a  rms  radius  of  about  1  cm. 

The  experimental  setup  which  has  been  used  to  observe  transverse 
oscillations^  in  the  ion  focusing  regime  is  shown  in  Figure  1.  The  drift  region 
is  constructed  of  15  cm  diameter  stainless  steel  tubes  with  a  total  length  of  1 
meter.  A  passively  integrated  Rogowski  coil  is  used  to  measure  the  net  current 
30  cm  downstream  of  the  diode.  A  magnetic  probe  array^  is  placed  adjacent  to 
the  Rogowski  coil.  Together  these  signals  will  give  information  about  the 
current  centroid  displacement  from  the  axis  as  a  function  of  time.  Further 
details  of  the  experiment  may  be  found  in  Reference  5.  In  order  to  perform 
open  shutter  photography,  approximately  60  cm  of  the  drift  tube  can  be  replaced 
by  a  Pyrex  tube  with  a  stainless  steel  screen  insert. 

The  experiment  is  performed  with  10,  20,  30,  40,  80,  or  160  mTorr  of  Argon 
or  Nitrogen  as  filling  gas.  In  all  cases,  transverse  displacements  are  detected 
from  the  probe  array  during  the  beam  pulse.  In  many  instances,  the  instability 
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is  so  violent  that  the  beam  current  centroid  is  displaced  a  large  fraction  of 
the  wall  radius  of  the  drift  tube.  In  several  cases,  at  the  lower  pressures,  it 
is  possible  to  identify  a  frequency  associated  with  the  oscillation.  The  signal 
at  higher  pressures  is  often  just  a  centroid  displacement  to  one  side.  Delay 
time  from  the  onset  of  the  current  pulse  at  the  probe  position  to  the  onset  of 
the  instability  was  also  observed  in  all  cases.  Figure  2  shows  traces  of  the 
beam  current  pulse  and  the  centroid  displacement  waveform  for  a  low  and  a  high 
pressure  discharges. 

In  order  to  relate  the  experimental  observations  to  the  theoretical 
calculations,  we  note  that  the  axial  wavelength,  xz,  of  the  perturbation  is 
given  by 


kz  =  2tt/xz  =  tto/IOO  cm  *,  n  =  1,2,3,  ...  (1.3) 

since  the  eigenfunctions  must  vanish  at  both  conducting  endplates  of  the 
1  meter  long  drift  tube.  Making  use  of  the  relevant  experimental 
parameters  (Ib  =  4  kA,  R  =  1  cm,  y  ~  2.4,  T^  =  35  keV)  and  Equation  (11),  we 
obtain  kz  >  0.254  cm"*,  which  indicates  xz  =  25  cm.  This  axial  wavelength  is 
evidenced  in  the  open  shutter  photograph  shown  in  Figure  3,  taken  with  the  Pyrex 
tube  in  place. 

The  r  ci llation  frequencies  for  very  low  pressure  cases  can  now  be  found, 


with  x  =  25  cm,  from  Equation  (9)  to  be  17.6  MHz  for  Argon  and  25  MHz  for 


atomic  Nitrogen.  These  values  are  in  good  agreement  with  experimental  data  at 
10  mTorr,  which  are  17  ±  3  Mhz  and  27  ±  1  MHz  for  Argon  and  Nitrogen, 
respect i vely. ^ 

For  both  Argon  and  Nitrogen,  the  characteristic  ionization  time  tQ  is  given 
approximately  by  1200/P(mTorr )  nanoseconds .  ^  Equation  (12)  can  then  be  used  to 
estimate  the  delay  time  for  the  onset  of  the  instability,  which  is 

*<!  5  PfmTorr )  nanoseconds.  (14) 

As  can  be  seen  from  Figure  4,  Equation  (14)  well  characterizes  the  behavior  of 
the  experimentally  determined  onset  time  for  Argon  (a  similar  situation  occurs 
for  Nitrogen). 

In  conclusion,  ion  hose  instability  of  an  intense  relativistic  electron 
beam  propagation  in  a  beam-induced  plasma  channel  has  been  investigated. 
Ex[_,imental  observations  agree  well  with  the  theoretically  calculated  values  in 
terms  of  the  axial  wavelength,  oscillation  frequency,  and  the  onset  time  of  the 
ion  hose  instability. 

This  work  was  supported  by  the  Independent  Research  Fund  at  the  Naval 
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ROGOWSKI  COIL 


50  ns/DIV 


OSCILLOSCOPE  TRACES  OF  THE  PASSIVELY  INTEGRATED  ROGOWSKI 
COIL  AND  THE  DISPLACEMENT  PROBE.  THE  ROGOWSKI  TRACE  (A) 

IS  1.2  kA/DIV.  THE  DISPLACEMENT  PROBE  TRACE  (B)  IS  0.94/KkA) 
cm/DIV  FOR  SMALL  DISPLACEMENTS  MUCH  LESS  THAN  THE  DRIFT 
TUBE  RADIUS. 

FIGURE  2.  CURRENT  AND  DISPLACEMENT  TRACE 
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PLOT  OF  THE  INSTABILITY  ONSET  TIME  AS  A  FUNCTION  OF  PRESSURE 
FOR  ARGON.  THE  ERROR  BARS  ABOUT  THE  EXPERIMENTAL  POINTS  ARE 
SAMPLE  STANDARD  DEVIATIONS.  THE  SOLID  LINE  REPRESENTS  THE 
THEORETICAL  PREDICTION  DERIVED  IN  THIS  PAPER. 

FIGURE  4.  INSTABILITY  ONSET  TIME  VS  PRESSURE 
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